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On the Symmetry of Low-Field Ordered Phase of PrFe 4 P 12 : 31 P NMR 
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We have performed 31 P nuclear magnetic resonance (NMR) experiments on the filled skut- 
terudite compound PrFe4Pi2 to investigate the exotic ordered phase at low temperatures and 
low magnetic fields. Analysis of the NMR line splitting in the ordered phase indicates that to- 
tally symmetric (Pi-type) staggered magnetic multipoles are induced by magnetic fields. This 
is incompatible with any type of quadrupole order at zero field. We conclude that the ordered 
phase has broken translational symmetry with the wave vector Q = (1, 0, 0) but the Th point 
symmetry at the Pr sites is not broken by the non-magnetic Ti-type order parameter. 
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Intermetallic compounds with the filled-skutteruditc 
structure RT4X12 (R — rare earth, T — transition metal, 
X = pnictogen) has attracted strong recent attention 
because of a variety of intriguing phenomena in a com- 
mon crystal structure such as metal-insulator transi- 
tion, 1 -' multipole ordering, 2 ' exotic superconductivity, 3 -' 
and anomalous phonons. 4 -' Among them, PrFe4Pi2 shows 
a peculiar phase transition at Ta=6.5 K. 5 ) The low tem- 
perature phase has no spontaneous magnetic moment at 
zero field 6 -' and is suppressed by a magnetic field of 4- 
7 T dependent on the field directions, 5 ' resulting in a 
heavy-fermion state with a large cyclotron effective mass 
m* = 81 mo- 7 ) The low temperature phase has a struc- 
tural modulation with the wave vector Q = (1,0,0), 
indicating loss of (|, I, |) translation. 8 ) Spatial order- 
ing of distinct electronic states of two Pr 3+ ions in the 
bec unit cell was also observed by resonant X-ray scat- 
tering, 9 ' and the field-induced staggered magnetization 
was observed by neutron scattering. 10 ) Although these 
experiments and the elastic measurements 11 -' suggest an 
antiferro-quadrupole order likely to be of I^-type, direct 
identification of the order parameter has not been made 
yet. 

In this letter, we report results of nuclear magnetic res- 
onance (NMR) experiments on 31 P nuclei (spin 1/2) in a 
single crystal of PrFe4Pi2, focusing on the ordered phase. 
We observed field-induced splitting of the 31 P NMR lines 
below Ta- The results for various field directions revealed 
that the splitting is due to totally-symmetric staggered 
magnetic multipoles (octupole as well as dipole) belong- 
ing to the Ti representation. We conclude that the order 
parameter at zero field has Ti symmetry, excluding any 
type of quadrupole order. 

In the filled skutterudite structure (the space group 
Jm3, T^), Pr atoms are surrounded by an icosahcdral 
cage of P atoms and form a body-centered cubic lattice as 
shown in Fig. I. 6 ' 12 ) Iron atoms sit at the middle between 
the corner and the body-centered Pr atoms. Although all 
the P sites are crystallographically equivalent, they have 



different NMR frequencies in magnetic field because of 
anisotropic hyperfine interaction. For later discussion, we 
define six types of P sites giving distinct NMR frequen- 
cies for general field directions as P1(0, u, v), P2(0, u, v), 
P3(u, 0, u), P4(v, 0, u), P5(u, v, 0) and P6(u, v, 0), where 
u and v are the asymmetry parameters (Fig. I). 13 ) 

If the direction of the external field H is invariant un- 
der a symmetry operation which transforms one P site 
to another, those two P sites must have the same NMR 
frequency. For example, PI, P3 and P5 sites (P2, P4, and 
P6 sites) transform each other by 120° rotation around 
[111], therefore, they should give a single NMR line for 
H || [111]. Likewise, for H || [001], there should be three 
lines generated by three pairs of sites, (PI, P2), (P3, P4), 
and (P5, P6). For H [| [110], P5 and P6 become inequiv- 
alent and we expect four lines. Figure 2(a) shows the 
NMR spectra at T=50 K above Xa- We observed three, 
four, and two lines for H parallel to the [001], [110], and 
[111] directions, respectively. This agrees with the conse- 
quence of crystal symmetry. Detailed angle dependence 
of the resonance frequencies (not shown) leads to the site 



[001] 




Fig. 1. (Color online) Schematic view of the crystal structure of 
PrFe4Pi2. The cage surrounding the body-centered Pr atom is 
omitted for clarity. 
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Fig. 2. (a) 31 P NMR spectra at 50 K obtained by Fourier trans- 
forming the spin-echo signals. The numbers indicate the P sites 
(PI— P6) to which each peak is assigned. The dotted line shows 
the reference frequency at zero shift, (b) 31 P NMR spectra at 
6.0 K (above T A ) and at 4.6 K (below T A ) for the field of 1.6 T 
along [111]. NMR lines from PI, P3 and P5 sites (P2, P4 and 
P6 sites) are indicated by solid (open) circles. 
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Fig. 4. (Color online) Field dependence of the splitting of the 
hyperfine field at 2 K for the field along (a) [001], (b) [110], and 
(c) [111]. The boundary between the ordered and uniform phases 
are shown by the dotted lines. 



A remarkable result is that AH extrapolates to zero at 
H=0. Such a feature has been also observed in poly- 
crystals by Ishida et al. 14 ) The field dependences of AH 
are different for different sites. This is most evident for 
H || [001], where AH for P3 and P4 increases mono- 
tonically with increasing field while AH for other sites 
shows a maximum. The field variation of AH is smooth 
in all cases without any jump or kink, indicating ab- 
sence of additional phase transitions up to the boundary 
to the high-field heavy fermion (HF) state. A first-order 



Fig. 3. Temperature dependence of the NMR frequencies near 
T A =6.27 K for the field of 0.4 T along [001]. 



assignment shown in Fig. 2(a). The small splitting of the 
two high frequency lines for H || [001] is probably due to 
nuclear spin-spin coupling. 

On crossing Ta from above, each line splits into a pair 
of lines as shown in Figs. 2(b) and 3. The splitting grows 
continuously near Ta at low fields (0.4 T) as shown in 
Fig. 3. At higher fields above about 1.5 T, however, the 
splitting develops discontinuously and there is a narrow 
temperature range in the vicinity of Ta where the spec- 
trum has both split and unsplit lines. This indicates that 
the second order phase transition at low fields changes 
to first order at higher fields, consistent with the results 
of specific-heat measurements. 5 -' 

We define the splitting of the hyperfine field AH as the 
frequency interval of the splitting divided by the nuclear 
gyromagnetic ratio 7 at. The field dependence of AH is 
plotted in Fig. 4 for various field directions at T=2 K. 



Table I. Symmetries of multipoles up to hexadecapoles in the 
crystal field. The + and — signs show the parity under 
time reversal. The multipoles are defined in terms of the dipole 
J as follows. 16 - 17 ) Quadrupoles: O!? = §(3J 2 — J 2 ), 

^(J 2 - J 2 ), 0(. v = ^J~~7^, Octupoles: T xyz = 
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subscripts. Since (T^,T^,T z a ) and (Tg,T y , 
and (Hx , Hy , H^)) have the same symmetry T7 
group, they can be mixed. 
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Table II. Symmetries of multipoles in the crystal field in magnetic fields. The numbers in parentheses in the fifth column indicate the 
equivalent P sites and the last column shows the number of NMR lines for a single P12 cage in the presence of the magnetic multipolcs. 
They should be multiplied by two to get the total number of lines in the antiferro-multipole ordered phase with Q = (1, 0, 0). 
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transition to the HF phase is apparent from the sudden 
vanishing of AH at the phase boundary. 

Detailed examination of angle dependence of the spec- 
tra below Ta reported in ref. 15 revealed that all of the 
NMR lines above T\ split always into two lines and AH 
never vanishes for any field direction. Obviously, the dou- 
bling of the NMR lines should be ascribed to the loss of 
(|, 5 , 1 ) translation and the distinct electronic states of 
the two Pr ions at the corner (PrI) and the body center 
(PrII) of the original bec lattice. 8 ' 9 ) This should divide 
each of the P1-P6 sites into two sublattices P1(I)-P6(I) 
and P1(II)-P6(II). The former (latter) belongs to the P 12 
cage surrounding the PrI (PrII) sites. We should stress 
that there is no additional line splitting, i.e., the number 
of NMR lines from each of these cages is exactly the same 
as the number of lines above Ta- This means that the Th 
point symmetry at both Pr sites is preserved below Ta, 
which is not compatible with any type of quadrupolc 
ordering. In the following, we provide more precise argu- 
ments on these points. 

The hyperfine field at P nuclei is determined by the 
spin density distribution of Pr-4/ electrons through the 
dipolar and the transferred hyperfine interactions. The 
very local nature of the latter interaction causes nuclei 
to couple not only to dipole moments but to octupolcs 
and higher order magnetic multipoles. This is because 
the local spin density near a P nucleus can be nonzero 
when a Pr ion has finite expectation value of a high order 
magnetic multipole, even if the dipole moment, or the 
spatial integration of the spin density, is zero. Such a 
possibility has been first recognized by Sakai et al. ls ^ in 
their analysis of the NMR data on CeBg. 19 ^ 

Magnetic and nonmagnetic multipoles up to hexade- 
capoles are presented in Table I as the basis of irreducible 
representations of the Th group. In magnetic fields, they 
are decomposed into sets of smaller number of basis for 
reduced symmetries as shown in Table II. The NMR line 
splitting AH is due to some staggered magnetic multi- 
poles with Q = (1, 0, 0). The vanishing AH at zero field 
indicates that these staggered multipoles also disappear 
at zero field. Therefore, the order parameter (OP) at zero 
field must be a nonmagnetic staggered multipole, which 
is even with respect to the time reversal and does not cou- 
ple to nuclear magnetic moments. As discussed by Shiina 



et al., 16 ' however, OP at zero field and the field-induced 
multipoles must belong to the same irreducible represen- 
tation of the symmetry group reduced by the magnetic 
field. Thus identification of the field-induced magnetic 
multipoles allows us to determine the symmetry of OP 
at zero field. 

Using the invariant form of the hyperfine coupling at 
the P sites in the filled skutterrudite structure, 15, 20 - ) the 
difference of the hyperfine fields at Pn(I) and Pn(II) 
(n=l-6), AH(ra) = (AH { x n) ,AH y 
n=l and 2 as 



Aiji™-*), is given for 
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C21 Jy + c 23 T° ± (c 22 r z + c 24 T°), 
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where T£ yz7 T|, and J| are the staggered octupole and 
dipole moments and is the hyperfine coupling con- 
stants. For terms with ±, the + (— ) sign should be taken 
for PI (P2). Expressions for P3 and P4 (P5 and P6) are 
obtained by applying once (twice) simultaneous cyclic 
permutations i-ti/->z^3;in the subscripts of the 
multipoles and AH X — > AH y — > AH Z — » AH X . Since the 
external field is much larger than the hyperfine field, AH 
is equal to AH • h, where h is the unit vector along the 
field direction. From eq. (1), we can determine the equiv- 
alent sites and the number of NMR lines in the presence 
of these field-induced magnetic multipoles as shown in 
the last two columns of Table II. For example, if a stag- 
gered component of J y or T y were induced by the field 



along [001], eq. (1) tells AH { Z 1} ^ AH { /> . Thus PI and 
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Fig. 5. Field-induced dipole J for H || [111], (a) When J is paral- 
lel to H. i.e. J c ^ and (J a , J b ) = in the notation of Table II, 
PI, P3 and P5 are all equivalent, (b) If J a or J b ^ 0, PI, P3 
and P5 become inequivalent due to loss of three-fold rotational 
symmetry about [111]. 
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Table III. Reduction of irreducible representations of the 
group in magnetic fields. 



zero field 
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P2 should give distinct NMR lines and there should be 
four lines from a single cage. 

As mentioned before, the number of NMR lines for a 
single P12 cage remains unchanged upon entering into 
the ordered phase: three, four and two for H || [001], 
[110] and [111], respectively. From Table II, we can con- 
clude that the field-induced magnetic multipoles must 
have the Ti symmetry for all field directions. In other 
words, only the totally-symmetric multipoles can avoid 
additional line splitting. For dipoles, this means that the 
induced dipole moment is always parallel to the exter- 
nal field, which is consistent with the neutron scattering 
measurements. 10 - ) This can be intuitively understood as 
illustrated in Fig. 5 for the case H || [111]. 

Thus the OP must be totally symmetric in magnetic 
fields. Inspecting Table II, we find there is no such 
quadrupole. For example, retype quadrupoles O xy , O yz 
and O zx induce magnetic multipoles which do not be- 
long to Ti for H || [001]. Quadrupoles 0% and 0\ of r 23 
type induce dipole moments perpendicular to the field 
for H || [111]. Both cases lead to more NMR lines than 
experimentally observed. Therefore, there should not be 
antifcrro-quadrupole order at zero field. 

One notices at this point that the OP at zero field 
should be also totally symmetric. We can see this for- 
mally by studying reduction of irreducible representa- 
tions of the T h group under the [001], [110] and [111] 
fields as shown in Table III. It is apparent that only the 
Ti representation at zero field has always a Ti compo- 
nent in magnetic fields irrespective of the field direction. 

At present, we cannot identify the detailed form of the 
OP because the number of NMR lines depends only on 
the symmetry of the OP in magnetic fields. The OP may 
thus include various multipoles of different ranks in gen- 
eral, as far as they are nonmagnetic and have the Ti sym- 
metry. Such a multipole order can be caused, for exam- 
ple, by alternate breathing of the icosahedral cages which 
preserves local Th symmetry at the Pr sites. This type 
of lattice distortion has been observed in the insulating 
phase of PrRu4Pi2, 21 - ) where the antifcrro order of the 
hexadecapole H° oc + Jy + has been proposed. 22 ) 
In fact, this is the only totally-symmetric nonmagnetic 
multipole among those listed in Table I. Quantitative 
analyses of the splitting AH as a function of direction 
and magnitude of the field will give further information 
about the OP such as the rank of a dominant multipole. 
Some phenomenological approaches are in progress, 15 ' 20 ) 
although a microscopic theory is needed to uncover the 
mechanism of this interesting phase transition. 



In conclusion, we have presented 31 P NMR data in a 
single crystal of PrFe4Pi2, confirming antiferro order of 
nonmagnetic multipoles at low fields. There exist field- 
induced staggered magnetic multipoles, which is compat- 
ible only with the totally-symmetric order parameter at 
zero field. The Th point symmetry at the Pr sites is thus 
preserved in the ordered phase, excluding any type of 
quadrupole order. 
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